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Abstract 

We report on an extensive numerical investigation of the Kardar-Parisi-Zhang 
equation describing non-equihbrium interfaces. Attention is paid to the de- 
pendence of the growth exponent fi on the details of the distribution of noise 
p{£,)- All distributions considered are delta-correlated in space and time, and 
have finite cumulants. We find that (3 becomes progressively more sensitive 
to details of the distribution with increasing dimensionality. We discuss the 
implications of these results for the universality hypothesis. 
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The Kardar-Parsis-Zhang (KPZ) equation although originally introduced as a 

model of non-equilibrium interface growth, has acquired a broader significance over the 
past decade as one of the simpler examples of a strong-coupling system. To this date there 
exists no systematic analytic procedure for determining the properties of the model for large 
values of the non-linear coupling. Additional interest in the model stems from its intimate 
mathematical relation to two other important systems, namely the noisy Burgers equation 
10], and the equilibrium properties of a directed polymer in a random medium (DPRM) [§. 

Most work on this problem to date has concentrated on determining the properties of the 
system (e.g. values of dynamic exponents) in the strong-coupling phase. In pursuit of this 
goal, the 'hypothesis of universality' (HOU) has been generally adopted; namely that details 
of the model representation of the KPZ equation should not affect universal quantities such 
as exponent values. In this Letter we report on an extensive numerical investigation of the 
KPZ equation which shows that exponents are sensitive to the precise form of the noise 
distribution, this sensitivity becoming extreme in higher dimensions. We regard our results 
as being strongly suggestive of a breakdown of universality in the KPZ equation, but of 
course, they cannot constitute a proof of this assertion. To set the scene for what is to come 
we devote one paragraph below to the general theoretical framework of the KPZ equation, 
and a further paragraph to the current state of knowledge. 

Denoting the interface profile by h{'x, t) (which is defined perpendicular to the d- 
dimensional substrate), the KPZ equation has the form 

dth = uV^h + X{VhY + ^ , (1) 

where ^(x, t) is a stochastic source, generally taken to be delta-correlated in space and time: 
= Tlx lit p(^(x, t)). The canonical choice for the distribution p is a gaussian, and the 
HOU is invoked to argue that any other choice (so long as it has finite cumulants) will 
lead to the same large-scale behavior. The prime goal is the computation of the exponents 
which characterize the evolution of fluctuations, along with the value of the upper critical 
dimension du, above which one expects the exponents to saturate at their mean-field values. 
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In the language of the KPZ equation, there are three exponents of immediate interest: 
Starting from a flat interface, the mean square fluctuations grow as W"^ = {h'^)c ~ t"^^, 
but will saturate for large times in a finite system of linear dimension L - the saturated 
value is expected to scale a.s W ^ L^. The fluctuations over the entire temporal regime 
may be conveniently described by the two-point correlation function C{r,t) = {{h{r,t) — 
h{0,t)y) ~ r'^^f{r^/t). The exponent z is often known as the dynamic exponent and gives 
the fundamental scaling between length and time. The three exponents may be reduced to 
one independent exponent by two scaling laws: zf3 = x ^^^d z + x = 2. The first law comes 
from the scaling behavior of C as r — > cxd, and the second is a consequence of invariance of 
the equation under an infinitesimal tilt transformation. It is also worth mentioning that the 
non-linear transformation w = exp(A/i/z/) produces the linear equation 

dtw = z/V^w + {X/iy)w C , (2) 

which corresponds to the equation for a {d + l)-dimensional directed polymer in a random 
potential — . 

Given the extreme difficulty of any analytic progress, most investigations of the KPZ 
equation have been numerically based. These investigations fall into two categories. The 
majority of numerical work consists in simulating microscopic models which under the HOU 
are assumed to have the same large-scale behavior as the KPZ equation. Such models in- 
clude Eden growth 0, polynuclear growth models 0, the restricted solid-on-solid (RSOS) 
model and its close relative, the hypercube stacking model The second numerical 



approach is that of direct integration of the equation itself [10,11 1. This entails some sub 



tleties of discretization which have only recently come to light [|T^ . However, there has been 
reasonable agreement between these various approaches as regards the values of exponents, 
at least in ci = 1 + 1 (where there exists an exact analytic result of z = 3/2 [0) and 
d = 2 + 1. In higher dimensions, numerical work becomes more difficult due to the smaller 
linear dimension of the systems (and the consequent early onset of finite size effects), but 
still there has been general agreement that although exponents may not be so precisely de- 



termined, there is no sign of a crossover to mean field values {z = 2) for dimensions up to 
d = 7 + 1 1^. An interesting analytic method which may be applied to the strong coupling 
regime is mode-coupling theory |M (although it is based on an ad hoc neglect of vertex 



renormalization) which seems to support du = 4. Also there has been recent work |T5| based 
on short-distance expansion techniques in the renormalization group (RG), which indicates 
the bound du < 4, whereas a mapping to directed percolation |]l6l suggests that du < 5. 



Our original motivation for this numerical study was to integrate a recently proposed 



discrete equation |12|, which was shown to correctly capture the strong-coupling behavior 
of the KPZ equation on a lattice. Only during the course of our work did we discover the 
noise sensitivity of the exponents, which is the focus of this Letter. Nevertheless, before 
presenting our results, it is useful to briefly describe this improved algorithm, and also to 
exhibit its close relation to the algorithm for zero temperature DPRM P], which has also 



been numerically studied in the past |17 



The key point concerning the discretization of the KPZ equation, is that one is only guar- 
anteed to capture the strong-coupling physics by discretizing the directed polymer equation 
(0), and constructing the discrete KPZ equation by the inverse transform hi = {v/X) In(wj). 
The time discretization necessarily introduces a two-stage process - i) pumping with the 
noise, and ii) relaxing with the deterministic part of the equation. Explicitly one has the 
discrete KPZ equation in the form 



hi{t + A) = hi{t) + {u /X)\n\l + {Au /a^) ^ [e 

^_ j nn i 



X{hj-hi)/u 



(3) 



where A and a are the grid scales for time and space respectively. Taking the strong-coupling 
limit A ^ oo allows one to write the relaxation stage of the above equation in the much 
simpler form 

hi{t + A) = max (hi, {L}) . (4) 

j nn j V / 

This corresponds exactly to the zero-temperature DPRM algorithm but written here in 



terms of the field hi rather than Wi. Note that there are no adjustable parameters in the 
strong-coupling algorithm, except for the functional form of the noise distribution p(i^) (since 
the grid scales, and the noise strength may be scaled away). 

We have implemented the above algorithm in dimensions d + 1 with d = 2,3,4, with 
a fiat surface {hi{0) = 0) as the initial configuration.. The only adjustable 'parameter' in 
our simulations has been the function p{C,)- This function was taken to be either gaussian 
PgiO ^ 6xp(— ^^/2), or of the form Pa{0 ~ (o" — |^|)", — cr < ^ < cr (where cr is adjusted to 
maintain unit variance for each choice of a) . On varying a the distribution pa interpolates 
through the forms: bimodal (a \ —1), top-hat {a = 0), triangular {a = 1), and finally 
distributions with very rapidly vanishing tails as a — * oo. Note that all distributions are 
strongly localized and have an infinite set of finite cumulants. [They are thus distinct from 
distributions with power-law tails, which are known to change the exponent values [Q.] 

Our simulations are performed on lattices of size L'^, with averaging over samples. 
The simulations were of the size (L = 2048, = 8) for = 2, (L = 200, = 8) for d = 3, 
and {L = 60, iV = 8) for = 4. Although these simulations are of the largest scale possible 
within our resources (a cluster of DEC alpha workstations), we are aware that specialist 
computational groups can improve on the precision of our results. However, the qualitative 
features of interest in this Letter are convincingly clear from our simulations. 

The interface width W is plotted as a function of time in Figs.l, 2 and 3, for dimensions 
d = 2,3 and 4 respectively. As mentioned earlier, this quantity is expected to grow as 
W ^ (so long as the dynamic length scale is much less than L). It is clear that there 
is a dependence of [3 on the value of a, this dependence becoming progressively stronger in 
higher dimensions. (This is the reason for not presenting data in dimension (1+1), for which 
all distributions yield a value of (3 close to the exact result of 1/3 [|1^]). The measured values 
of (3 are presented in Table 1. An important point to make is that in all dimensions, our 
results are in agreement with the consensus of exponent values in the literature PJ9|,p!7[1, if 
we use the gaussian noise Pg{C)- However, the values of (3 drops smoothly as the parameter 
a is decreased. 
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Naturally one may attempt to interpret these results in terms of a temporal crossover. 
Given the straightness of the curves, such a crossover would be exceptionally slow. We have 
checked for the presence of this phenomenon by trying to fit the data with the generic form 
W'^ = At^^ + Bt^"' (which implicitly includes the popular fitting Ansatz 7 = used by 
previous groups P,P^). The best fits correspond to values of [3 close to the value one would 
obtain by the simpler fit ~ , with the 'correction to scaling' exponent 7 taking a value 
of approximately (3/2. This indicates that the curves show no sign of crossover, as they have 
a clearly dominant power-law form, with modest corrections. The quoted errors in (3 (shown 
in Table 1) are estimated from the range of (3 with which one may obtain an acceptable fit 
using the above fitting Ansatz. 

To enable a check of our results, we have concentrated on the value a = 1/2 m. d = ?> for 
which we averaged over = 64 samples of size L = 180. In this case, the data is of a good 
enough quality to measure a running value of the growth exponent Pes(t) = d\Ti{W) / dhi{t) 
||17|| . The quality of the data obtained for the value a = —1/2 (with A^ = 8 samples) is 
also sufficiently good to allow this measurement. The running exponents (3es{t) for a = 1/2 
and a = —1/2, plotted against t^^ (which magnifies potential systematic deviations from a 
simple power law), are shown in the insets of Figs. 4 and 5 respectively; with the previously 
measured values of (3 used in each case. There is no sign of any asymptotic deviation away 
from these values of (3 = 0.14 and (3 = 0.05. 

Another independent check was made for these two values of a, by studying the variation 
of the saturated (or steady-state) width Wss, as a function of the system size. On scaling 
grounds one expects Wss ~ L^, with x = 2/5/(1 -|- /5), where we have assumed the exponent 
relations zP = x -2 + X = 2. In Figs. 4 and 5, we plot Wss against L for a = 1/2 
and a = — 1/2 respectively. One observes that there is more than a decade of clean scaling 
in each case, with fitted values x = 0-24 and x = 0.08 respectively, which are consistent 
(within the stated errors) with the previously measured values of (3. On the grounds of this 
numerical work we are led to the statement that the growth exponent (3 is non-universal 
with respect to changes in the noise distribution p{$,). In the remainder of the paper we 
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present a brief discussion of the robustness and physical imphcations of this claim. 

It is possible to gain information about the exponents from physical quantities other 
than W. As an example, one may measure the two-point correlator C{r,t) (or its Fourier 
transform, the structure factor) and attempt to collapse the functions using dynamical 
scaling. We have been able to produce tolerably good data collapse for C{r,t) in all cases, 
yielding values of the exponents within the errors of those listed in Table 1. However, in 
our opinion the measurement of exponents from data collapse is less reliable than direct 
measurements of W due to the difficulties of including corrections to scaling in a systematic 
fashion. 

As mentioned above, one may try to interpret these results in terms of slow temporal 
crossover. Our analysis indicates that strong crossover effects are absent from the data, since 
we are able to measure reasonably precise values of the exponent f3 along with its 'correction 
to scaling exponent' 7. However, one can not rule out the possibility that one has simply 
failed to reach the 'true' asymptotic regime (AR), and that one is measuring some transient 
scaling regime. This introduces the question of how one should empirically define the AR. 
For our purposes we have used the standard working definition that in the AR one observes 
clean power-law behavior of the quantities of interest (in our case the interface width W), 
and that the correlation length in the system is large compared to the lattice scale. Both 



criteria have been met in our simulations. A further criterion which one may invoke [T^ is 
that the interface width itself must be much greater than the effective lattice spacing in the 
growth direction (here set to 0(1) by the noise variance being normalized to unity.) This 
criterion is not satisfied in our simulations for c? = 4, or for very small a, for the simple 
reason that if the interface fiuctuations grow very slowly with time (meaning f3 is small) 
then one will never achieve W ^ 1 on observable time-scales, even though the correlation 
length of the system is large (since z is generally close to 2). There can be no definitive 
answer to the question of whether a given simulation has reached the AR, although we would 
encourage specialist computational groups to improve on our time scales and system sizes 
in order to shed more light on this question. 



As a final comment on the simulations, it is possible that these results are a consequence 
of working at the strong-coupling limit (A oo). An example is known in the directed 
polymer literature [|1^] of exponents dependent on the noise distribution at zero temperature 
- for the case of perfectly correlated disorder in the longitudinal direction. The reasons for 
this dependency are physically clear and may be illustrated by a simple Flory argument. 
However, it is doubtful whether a similar effect is relevant to our simulations, in which the 
noise is delta-correlated in time. 

If the HOU is false for the KPZ equation, it is crucial to trace the underlying physical 
reasons for this. The most fragile aspect of KPZ physics is the role of the microscopic cut- 
off. For instance, the fact that the naive (but canonical) discretization fails to retain the 
continuum physics |jl2[ gives one real cause for concern. It may be that the lattice scale is 



always relevant to the large-distance scaling of the interface, which would then give scope 
for non-universal features, such as those seen in the present work. A similar example of non- 
standard scaling is known from controlled calculations on the deterministic Burgers problem 
|2§. 

In this paper we have given convincing numerical evidence that the KPZ growth exponent 
P is strongly sensitive to certain details of the noise distribution (here characterized by the 
parameter a). It is important to understand whether the results presented here are really 
the hallmark of a breakdown of the HOU for the KPZ equation, or whether there exist 
(exponentially) long crossover scales deep within the strong-coupling phase itself. Whichever 
is the case, we believe that these results are indicative of unexpected and interesting new 
physics in the KPZ problem. 

The authors are grateful to Alan Bray, Joachim Krug, Amos Maritan and Michael Moore 
for useful conversations. TJN acknowledges financial support from the Engineering and 
Physical Sciences Research Council. 
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List of figure captions 



Fig. 1: Interface width versus t in dimension 2 + 1. The upper curve corresponds to Pg 
and the lower to Pa with a = 0. The straight lines are fitted with values of /3 given in Table 
1. 

Fig. 2: Interface width W'^ versus t in dimension 3 + 1. The uppermost curve corresponds 
to Pg and the lower curves to Pa with a = 2, 1, 1/2, and —1/2 in descending order. 

Fig. 3: Interface width versus t in dimension 4 + 1. The uppermost curve corresponds 
to Pg and the lower curves to Pa with a — 1, 0, —1/2 in descending order. 

Fig. 4: Steady-state interface width Wgg versus system size L in dimension 3 + 1, using 
distribution p^ with a — 1/2. The fitted line has a slope of 2% = 0.48. The inset shows the 
running exponent Pesit) versus with /3 = 0.14. 

Fig. 5: Steady-state interface width Wgg versus system size L in dimension 3 + 1, using 
distribution pa with a — —1/2. The fitted fine has a slope of 2% = 0.16. The inset shows 
the running exponent (3ea{t) versus with /3 = 0.05. 

Table. 1: Values of the measured growth exponent /3 as a function of dimension d, and noise 
distribution p(^). 
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Figure 4, Newman and Swift, 'Non-universal exponents in interface growth' 
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Figure 5, Newman and Swift, 'Non-universal exponents in interface growth' 



